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Instanton-dyons, also known as instanton-monopoles or instanton-quarks, are topological con¬ 
stituents of the instantons at nonzero temperature and holonomy. We perform numerical simula¬ 
tions of the ensemble of interacting dyons for the SU(2) pure gauge theory. Unlike previous studies, 
we focus on back reaction on the holonomy and the issue of confinement. We calculate the free 
energy as a function of the holonomy and the dyon densities, using standard Metropolis Monte 
Carlo and integration over parameter methods. We observe that as the temperature decreases and 
the dyon density grows, its minimum indeed moves from small holonomy to the value corresponding 
to confinement. We then report various parameters of the self-consistent ensembles as a function of 
temperature, and investigate the role of inter-particle correlations. 


I. INTRODUCTION 

QCD description of strongly interacting matter at fi¬ 
nite temperature T has originated from 1970’s. At first, 
its high temperature phase - known as Quark-Gluon 
Plasma, QGP - has been studied using weak coupling 
methods, see e.g. reviews mm- 

The interest then switched to non-perturbative phe¬ 
nomena, related with the topological solitons of various 
dimensionality and two basic non-perturbative phenom¬ 
ena: confinement and chiral symmetry breaking. Instan¬ 
tons [3], the Euclidean 4-dimensional topological solitons, 
have at high T the sizes p ^ 1/T and appear with the 
probability 

"^instantons ^ €Xp[ StT fg (T)] ^ 

where the power is the one loop beta function coefficient, 
b = llNc/3 for SU{Nc) gauge theory. So, at high T 
the density is small and the topology is not important. 
Conversely, as T decreases, the instanton density grows 
rapidly, till they become an important ingredient of the 
gauge fields in the QCD vacuum. Chiral anomalies in¬ 
duce existence of the fermionic zero modes of instantons, 
which generate the so called’t Hooft effective interaction 
of 2Nf fermions, which explicitly violates the 17^(1) chi¬ 
ral symmetry. Furthermore, collectivization of the zero 
modes create the so called Zero Mode Zone of quasi-zero 
eigenstates, which break spontaneously the SU{Nf) chi¬ 
ral symmetry. Although those states includes only tiny 
subset of all fermionic states in lattice numerical simula¬ 
tions, they are the key elements of the chiral symmetry 
breaking and the hadronic spectroscopy. The so called 
Interacting Instanton Liquid Model (IILM) has been de¬ 
veloped, including ’t Hooft interaction to all orders, for 
a review see |3]. 

As the temperature decreases from the high-T regime, 
another important phenomenon is appearance of non¬ 
trivial expectation value of the Polyakov line. For the 
SU(2) gauge theory we will be discussing in this work, 
it is related to the so called holonomy parameter by 
(P) = cos(tti>) (for explicit notations see Appendix [A| . 


While at high T it vanishes jz —>■ 0, at temperatures at 
and below the critical value it reaches the so called 
“confining value” v = 1/2 at which the Polyakov line 
vanishes. This leads to switching out quark and gluon de¬ 
grees of freedom, and transition from QGP to hadronic 
matter. Study of the instantons at nonzero holonomy 
has lead Lee,Lu,van Baal and Kraan u m to the so 
called KvBLL caloron solution, which revealed that at 
It ^ 0 the instantons get split into A7 (number of col¬ 
ors) (anti)dyons, (anti)self-dual 3d solitons with nonzero 
(Euclidean) electric and magnetic charges. (Details are in 
Appendix]^. Because of long-range forces between these 
dyons, we will thus refer to instanton-dyon ensemble as 
the “dyonic plasma”. 

Unlike instantons, the instanton-dyons interact di¬ 
rectly with the holonomy. Diakonov [7j suggested that 
back reaction of the dyon free energy on holonomy is 
responsible for confinement phase transition but was un¬ 
able to show it. Poppitz, Schaefer and Unsal [5] had 
shown that instanon-dyon confinement does occur in a 
very specific “controlled setting”, a supersymmetric the¬ 
ory compactified on C with a small spatial cir¬ 
cle and periodic fermions. The smallness of the circle, 
like high T, makes the coupling weak. The periodic 
fermions preserve supersymmetry and cancel the pertur¬ 
bative holonomy potential Vgpf(j^), which allows con¬ 
finement to be induced even by an exponentially small 
density of the dyons. These authors have been able to 
trace the crucial effect to the repulsive dyon-antidyon in¬ 
teraction inside the dyon-antidyon pairs, which they call 
“bions”. 

A phenomenological model showing that repulsive in¬ 
teraction between them, modeled by an excluded volume, 
has been proposed for QCD-like theories by Shuryak and 
Sulejmanpasic [9], which reached qualitative description 
of the deconfinement phase transition and other proper¬ 
ties of the thermal SU{2) pure gauge system above Tc, 
in qualitative agreement with available lattice data. We 
will discuss similar model in section |llj before we embark 
on numerical simulations. 

Although the interaction between the instanton-dyons 
have been studied for a long time, the leading-order ef- 
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feet - classical dyon-antidyon interaction has been miss¬ 
ing. The corresponding studies, deriving the so called 
“streamline” set of configurations via the gradient flow 
method, has been done in our previous work |10j . This 
classical interaction turns out to be weak in relative terms 
SS S but, still being classical SS ^ 1/g^, rather large 
numerically SS = 0(1) to induce significant correlations 
in the ensemble. 

Liu, Shuryak and Zahed m have recently shown that 
one can incorporate this classical interaction by the mean 
field techniques, but only if the ensemble is dense enough 
to generate sufficient screening. In terms of the temper¬ 
ature, this treatment applies only for T < T^. 

The goal of our present work is to study the instanton- 
dyon ensemble by the direct Monte-Carlo simulation, 
without the mean field or any other approximations. As 
we will show, this will allow us to cover both a dilute and 
a dense regimes, and see in details the transformation of 
the holonomy potential which drives the deconfinement 
transition. In a way, this work is complementary to El, 
since our main focus is at temperatures T > Tc- 

Technically, the details of the setting to a large extent 
follow the first Monte-Carlo simulations of the instanton- 
dyon plasma by Faccioli and Shuryak HU. One major 
difference is the inclusion of the classical “streamline” 
interaction which were not known at the time of that 
work. The other is that that paper focused on the role 
of fermions and chiral symmetry breaking rather than 
confinement. (We expect to report on our next paper, 
with fermions, soon.) 

The paper is structured as follows: standard informa¬ 
tion about our notations, the holonomy potential and 
the instant on-dyons are delegated to sections of the Ap¬ 
pendix. We start in section |TT] by introducing a simple 
model which illustrates the main physics under discus¬ 
sion. Then in section III we explain the setting of the 
simulations, the interaction between the instanton-dyons 
and the moduli spaces which provides the measure in 
the partition function. In section |IV| we describe how 
we make the actual simulations and evaluate the free en¬ 
ergy. The back reaction of the ensemble on the holonomy 
potential is described in section |Vj which is followed by 
“self-consistency” study of the parameters in the section 
VE The physical results are summarized in the section 
VII| those include the holonomy potential and the screen¬ 


ing masses, as well as the densities of all types of dyons. 


II. AN EXCLUDED VOLUME MODEL 

To understand the main physics involved and the qual¬ 
itative behavior of the ensemble, including the confine¬ 
ment phase transition, we start by a discussion of a sim¬ 
plified model in which the only interaction is the repulsive 
core, making the volume occupied by each particle un¬ 
available to others. It is similar in spirit to that proposed 
by Shuryak and Sulejmanpasic |U, but is somewhat closer 
technically to the simulations to follow. 


We work in dimensionless units and we therefore define 
the 3-volume the density nt = ^, and the 

free energy density as = /. More information on 

units and notations can be found in Appendix. 

The effect of the excluded volume is accounted for in 
a very schematic way, by cutting off the partition func¬ 
tion when the amount of available volume vanishes. The 
volume of the M and L dyons scale by I/iz^ and Xjifi 
respectively, with P = 1 — v. We thus have the partition 
function in dimensionless units as a sum limited from 
above 


V3/{Vo)<M/u^+L/p^ 

E 

M,L 


exp 


~ 47r^ 2-2 

— Vs- u u 


X 


du — 

S = 


exp(-5'j^) 

Stt^ 


( 2 ) 


(3) 

(4) 


If the upper limit is ignored, and the volume goes to infin¬ 
ity, log{Z)/V^ —^iy‘^D'^ + 2{d^+dij), the perturbative 
homonomy potential plus the contribution of the nonin¬ 
teracting dyons. This explains that di, - the semiclassical 
dyon amplitude - is their density if all interactions are 
ignored. The parameter S is in fact the classical action 
of the caloron, or L + M system. The square comes from 
assuming the same amount of dyons and antidyons. 

If we are in the confining phase, v = v = 1/2, all 
dyons have the same sizes, and it is easy to introduce the 
excluded volume, for N dyons via 

—>■ t/j (Vs — Vexcluded)---{V3 — {N — l)Vexcluded^ 


However, since in general we have L, M dyons of differ¬ 
ent sizes, the analogous expression becomes cumbersome. 
Experimenting with those, we observe that similar results 
are obtained by simply cutting out the sum at the num¬ 
ber when there is no volume left, V 3 < Vq{M/ v^ + L/v^). 
where Vq is the excluded volume normalized for a dyon 
dX V = 1. 

Using Sterlings formula nl « 3 j 2 'im ((() for a large 
volume, we rewrite the sum as 


V 3 /(Vo)<M/iy 3 +L/p 3 

Z = ^ exp 

M,L 


2nM In 

d^e 

— 2rzi In 

dpc 

)1 


nM_ 


1 . '^L J 

J\ 


The free energy given by F{h') = —TlogZ depends on 
zz, located in the cutoff, in the dyon parameter dy, and in 
the last perturbative term. This last term, if dominant, 
would select z/ = 0 or z? = 0. (The value of S, the action 
for a dyon with v = I, is treated as an input parame¬ 
ter.) If one wants the two former (dyonic) exponents to 
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FIG. 1: (Color online). Plot of the free energy density / 
as function of holonomy v for A = 0.5 and Vb = 0.3 at y = 
4, 3.5, 3, bottom to top. It is seen how the maximum as a 
function of g goes further and further towards the conhning 
value of 1/2 as g goes up, and S and T go down. 


be dominant instead, and select some nontrivial v the 
corresponding densities should be large enough. In the 
simplest confining case ^ = P = 1/2 and dn = dp, one 
finds that the density has to be of the order of 1 /4 of the 
perturbative energy to overpower it. 

The expression ([^ is put into Mathematica and the 
maximum is found, for large enough volume, say V = 
900. One finds a sharp peak in N distribution, defining 
the density. Finding the maximum as a we vary we 
get f{v) = —logZ/Vs plotted in Fig. [^for g = 4, 3.5,3. 
At smaller g (larger S and higher T) the dyons are more 
suppressed and the free energy density / has a minimum 
at smaller ly. For increasing coupling g (decreasing S 
and T), the minimum shifts from zero, eventually to its 
confining value v = 1/2. 

Instead of showing z/ itself, one can also plot the aver¬ 
age Polyakov loop, depending on it as < P >= cos{'kv), 
versus S in Fig. The parameter S grows with T (see 
details in g. In this model it is seen how the Polyakov 
loop continuously goes to 0 at S' slightly smaller than 6. 
In order to get a better perspective on how this happens, 
we show the free energy versus the holonomy, as three 
curves in Fig[T] Note that in the critical case, the middle 
curve, the free energy f{iy) becomes very flat. 

As the repulsion (excluded volume) grows, the free en¬ 
ergy increases, and start coming into the positive do¬ 
main, see Fig|^ The / > 0 formally means that the 
vacuum without any dyons, with / = 0, is lower than 
with them: too repulsive dyons drive themselves out of 
existence. Obviously such situation is not physically pos¬ 
sible. 

We then show in Fig. at above the deconfinement 
transition the densities of different type {M and L) dyons 
is different. We will see similar plots as a result of numer¬ 
ical simulations below. Here let us only comment that di¬ 
rect evidences for um > in the deconfined phase have 
been found on the lattice, and are related to the issue 
of chiral symmetry breaking for fermions with different 
periodicity angle, see discussion in [15] . 



FIG. 2: (Color online). Polyakov loop P as a function of 
action parameter S for A = 0.5 and Vo = 0.3. 



FIG. 3: (Color online). Densities rii oi i = M or i = L 
dyons as a function of the action parameter S, for A = 0.5 
and Vq = 0.3. Note that the two densities are different at 
S>6. 



V 

FIG. 4: (Color online). The same as Fig[T] but for 

larger excluded volume Vo = 0.6 (and slightly different g = 
4,3.5,3.25). Since f — 0 corresponds to no dyons, a jump 
will be observed from confining holonomy v = 1/2 to trivial 
holonomy v = 0. 


III. THE INSTANTON-DYON INTERACTIONS 

Important new element is the inclusion of the leading 
order dyon-antidyon interaction recently studied in our 
previous paper m- We will use a slightly different pa- 
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rameterization that follows from the same data 

ASdd = - 1.632e-° ™“) 

9^ X 

X = 2Tn/rT 


( 6 ) 


for distances larger than x > A. 

It has been found that at interdyon distance x = 4 
the streamline terminates by rapid annihilation of the 
magnetic charges. If dyons are put at smaller distances, 
they repel till distance 4, before annihilation. Those con¬ 
figurations were not yet studied in detail, and thus our 
potential for x < 4 is a reasonable guess. We cut the po¬ 
tential off at a distance xq = r^T [2ttv) with a core which 
we describe by 




lyVo 


1 -I- exp [a{x — xo)] 


(7) 


where we scale the core also by v since we want the in¬ 
teraction to disappear at j/ = 0 . 

It is perhaps worth reminding why the attractive 
Coulomb-like interaction is there. Selfdual and antiself¬ 
dual objects have electric and magnetic forces canceling 
each other, but the effect still comes via A 4 and the non¬ 
linearity of the field strength tensor. It is well known 
from the analytic results for calorons that such inter¬ 
action can be described purely from the charges in ta¬ 
ble |l] such that the interaction is C = (6162 -I- mim 2 — 
2 ^ 1 / 12 ) Note that this interaction is repulsive for 

ML,LM channels, a result that has been checked by 
us (after submission of our previous paper [IH])- Note 
also that this expression also agrees with zero classical 
interaction between sectors that are completely self-dual 
(LM) or anti-self-dual. 

The volume element of the space of collective variables 
we will use is in the form of the so called Diakonov de¬ 
terminant 


y/g = detG 


( 8 ) 




k^i 


^ T\Xi_^rn ^k,r. 


(9) 


+2E 


T\Xi^rn ^k,p^r 


r) 


+2(5^7 


rTi\ 

1 ^ 2 ,m 


- 25, 


'm^n 




T\Xi^rn 


where Xi^m denote the position of the i’th dyon of type m. 
The Diakonov determinant is a combination of the metric 
between a M and L dyon, which is true at any distance, 
and the metric between dyons of same type, which is only 
true at large distance, but does have the correct behavior 
of being repulsive for small separations. We therefore 
introduce a cutoff on the separation r + cutoff^, 

such that for one pair of dyons of same type, the diagonal 
goes to 0 for iz = 0.5, instead of minus infinity. 


We use the same metric for the antidyons also. 

It is observed that in case the density of M and L 
dyons are different, then in the infinite volume limit, the 
sum will diverge. We therefore regularize all the terms 
with r — >■ where we work with the dimensionless 

Debye mass. 

The same exponentially damped effect due to the De¬ 
bye mass is also applied to the Coulomb potentials. With 
this the interaction is given by 

(^(eie 2 - 2/11/12)--l-mim 2 -^ 

9^ \ X xj 

X = 2'Ki'rT (10) 

for r larger than the core of size ro/(27rjz) for all com¬ 
binations except between dyons and their antidyon. For 
the dyon antidyon potential we have 

ASjjo = -2?^(-- 1.632e-° ™“)e-^^’'^ 

9 X 

X = 2'Ki'rT (II) 

We include the core for both dyon antidyon interaction, 
but also for dyon dyon interaction due to the lack of a 
repulsion, which otherwise destroys the simulation. We 
hope that such an interaction can be found due to loop 
calculations as it was done with instantons. 

= - -r rpf -^ ( 12 ) 

1 -I- exp [crJ (r — ro)(27riz)J 


IV. THE SETUP 

Like in |12j . instead of toroidal box with periodic 
boundary conditions in all coordinates, our simulations 
have been done on a sphere (in four dimensions), to 
simplify treatment of the long range Coulombic forces. 
In this pilot study we fix the total number of dyons to 
64. We opt instead for multiple runs, displaying depen¬ 
dencies of the free energy on all parameters, sacrificing 
somewhat statistical accuracy. We do not use large com¬ 
puters, relying on multiple cores of standard GPU’s of 
one standard computer. 

The radius of the sphere together with the ratio of M 
dyons to L dyons have been used to change their density. 

Iteration of the system is defined as a loop in which 
each dyon has had its position changed and the new 
action has then been accepted with the probability of 
exp(—A^) via Metropolis algorithm. The typical num¬ 
ber of iterations, for equilibration is 400 and after equi¬ 
libration 1600. 

In order to get the free energy we also use standard 
method. One can differentiate with respect to an aux¬ 
iliary parameter A introduced in front of the action and 
get 

= J Dx exp(—AS'(x)) (13) 

= T{S) 


g-FW/T 

dF 


(14) 
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Since the free energy at A = 0 is known analytically, one 
can integrate up to get the free energy at A = 1. When 
we do this we of course need to be careful about regions 
with a quick change in the action. 

For the calculation of det G it has been observed by 
Bruckmann et al [14] that it is only make sense if all 
eigenvalues are positive. It was observed El that, for 
randomly placed dyons this is typically not the case, un¬ 
less density is very low. In HH this issue has been dis¬ 
cussed further, with a conclusion that the Diakonov de¬ 
terminant can remain positive definite at higher densities 
needed, but only with certain correlations in the dyon lo¬ 
cations enforced. 

We have therefore used the Householder QR algorithm 
together with tridiagolization of the matrix G El to find 
the eigenvalues. We also redefine the potential as follows: 

If all eigenvalue are positive 


to dyons’s collective variables. 

^ — ^unchanged^changed 

The weight for one caloron {L M pair) was explicitly 
calculated in [I6j : at zero holonomy it agrees with the 
instanton result by ’t Hooft. Part of the answer is the 
factor coming from metric volume element ^fg in the 
space of L,M collective variables. Later Diakonov |7] 
combined this result with previously known answer for 
metric of two monopoles of the same kind (e.g. M, M 
pair) into an elegant expression for any number of L, M 
dyons now called Diakonov determinant det G. Taking 
the dilute limit ri 2 —> oo in both cases both formulas 
reduce to the same ri 2 dependence and one finds that 
the caloron weight from m needs to be divided by the 
factor {4:Tri'){4'KD) (see appendix 


Vd = - \og[Det{G)] : Vd < (15) 

Vd ^max • ^ ^max (16) 

and for one or more negative eigenvalues 

Vd = Vmax (17) 

This is done such that the excluded volume from the 
regions of negative eigenvalues are included, and at the 
same time to not create a region where the configuration 
is confined inside the region of negative eigenvalues. The 
possible great jump in energy from this definition means 
that we first integrate the free energy up to A = 0.1 
finely with 10 points, and then with 9 points integrate 
from A = 0.1 to A = 1. 


V. THE DYON BACK REACTION: 
HOLONOMY POTENTIAL 


In the absence of non-perturbative effects, there is only 
a perturbative interaction of thermal gluons with the 
holonomy, generating the perturbative Gross-Pisarski- 
Yaffe potential |5] (A2) which disfavor confinement. 
However, as reproduced countless of times in any finite 
T lattice gauge simulations during the last 3 decades, 
the peak of the holonomy distribution shifts to confine¬ 
ment at T « Tc- The corresponding effective potential 
V{v) has been numerically studies and parameterized, 
used in many models of finite-T QCD, e.g. the so called 
Polyakov-Nambu-Jona-Lasinio (PNJL) model. Now our 
task is to derive such potential, stemming from back re¬ 
action of the instanton-dyons. 

So, we add Vqpy to the dyon free energy obtained from 
our simulations and determine the total free energy of the 
system (obviously, assuming that there are no other rele¬ 
vant non-perturbative contributions). The dyon-induced 
partition function is further split into two factors: one 
containing all factors which depend on parameters un¬ 
changed in a simulation sets, and the second one related 


^unchanged — I ^2 J ^ 9 1/3 j/ 3 (19) 

Note that at the trivial holonomy i/ ^ 0 limit Z^^changed 
is ~ l/v. it is to be canceled by the diagonal part of the 
det(G) in the second part. 

Since we want to do the simulation for different amount 
of M and L dyons we need to split up the weight into a 
M part and L part and sum over all number of particles. 
We choose 


■^unchanged 


E 

Nm,Nl 



( 20 ) 


NmI ^ 


X 


1 




Nl 


2 


where we use that the amount of dyons and antidyons is 
the same. We simplify this as 




unchanged — 


Nm,Nl 

1 


/L 

! (l^*) 


Nm 


dj, = A ( ) e i/°3 V(47r) (21) 


Zchanged IS the interactions explained in section HI and 
thus also depends on the number of particles 




changed 


A/ = 


yml+NM) I D^^det{G)exp{-ADjgDix)) 

^^E>{^changed^ 
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where we have normalized such that Zchanged = 1 for no 
interaction. Combining with the unchanged part and the 
purturbative potential we get in the limit y —>■ oo 




H exp 
Nm,Nl 



47r^ 9 9 

d^e 

-V 3 

- 1 / — 2nM In 


\ 

L 3 

Um _ 


+ 2 nL In 


dpe 

riL 



(23) 


For 1^3 —>■ oo the partition function is completely domi¬ 
nated by the maximum of the exponent. Finding the free 
energy corresponds to finding the minimum of 


/ 


47r^ 

IT 


2 -2 
V V 


— 2 nM In 


d^e 

Um 


— 2 nL In 


dpe 

riL 


+ A/ 


(24) 


Note that as the dyon density increases, it changes 
its shape, producing a non-trivial minimum at v ^ 0 . 
Furthermore, at high density this minimum moves to iz = 
1 / 2 , the confining value. 

The densities of both kinds of dyons riL^riM are not 
in general equal: the model should be able to do this by 
adding compensating charge to the whole sphere. In our 
model this is done by including the Debye mass. 


VI. SELF CONSISTENCY 



V 

FIG. 5: (Color online). Free Energy density f as a function of 
V at S = &, Md = 2 and Nm = Nl = 16. The different curves 
corresponds to different densities. • n = 0.53, M n = 0.37, 
♦ n = 0.27, k n — 0.20, T n = 0.15, o n = 0.12. Not all 
densities are shown. 


changing the holonomy which is the Debye mass. We 
require that the used value for the Debye mass is the 
same as the one found from the derivative of F^ or atleast 
not more than 0.4 below the used value. 

The results shows that as the Debye mass goes to zero 
around the phase transition the only configuration that 
is consistent with this is that of equal M and L dyons. 


The partition function we simulate depends on several 
parameters, changed from one simulation set to another. 
Those include (i) the number of the dyons (ii) 

the radius of the sphere r; (iii) the action parameter 
S; (iv) the value of the holonomy i/, (v) the value of the 
Debye mass Mjj; (vi) the auxiliary factor A, which is then 
integrated over as explained in section |IV| 

In principle, the aim of our study is to obtain the de¬ 
pendence of the free energy on all of those parameters (i- 
v). While the practical cost of the simulations restricts 
the number of points one can study, we still had gen¬ 
erated more than hundred thousand runs and multiple 
plots. However, most of it neither can nor should be 
included in the paper. Since our physics goal is to un¬ 
derstand the back reaction of the dyon ensemble on the 
holonomy, we study the whole range of holonomies, from 
v = 0 to v = 1 / 2 , and only then locate its minimum. As 
for the Debye mass, we will find it from the potential and 
then show only the “selfconsistent” input set. 

What we actually need to describe at the end is not 
the free energy in the whole multi-dimensional space of 
all parameters, but the location of the free energy min¬ 
ima. The resulting set should be of co-dimension 1, since 
the original physical setting of the problem - the gauge 
theory at finite temperature - has only one input param¬ 
eter, T. 

2 ^^2 pi cy 

Using the definition of the Debye mass ^ ^ 2/7 = Mfy 
for fixed density we get the configurations response to 


VII. THE PHYSICAL RESULTS 

We now show only the result which fulfill the self- 
consistency requirement. Without fermions the results 
are symmetric in v ^ 1 — v and the results are therefore 
only for v < 1/2. We have included the Diakonov deter¬ 
minant, though its impact is not too great due to the not 
so small Debye mass which has been calculated using 3 
points. The results here are shown for a wall of 2/[2'kv) 
which was chosen in order to have a large enough density 
of dyons to overcome the purturbative potential, without 
completely making the perturbative potential irrelevant. 
We used A = 1.5 to obtain a phase shift around S = 6 . 
Action is related to temperature as explained in appendix 
[Aj This should of course be fitted to numerical data, but 
the present data on dyons does not have a high enough 
efficiency of detection to do this. The action goes up to 
S = 13, beyond this value the number of L dyons become 
too close to 1 , and we would need a higher total of dyons 
to proceed. 

The first thing to note about the results is that due 
to the repulsive Coulomb term between dyons and an- 
tidyons of different type, the free energy preferred to have 
a large Debye mass due to cutting off this repulsion. This 
meant that when the free energy spectrum as a function 
of holonomy for a fixed density becomes flat, the small 
Debye mass created a rise in energy. This resulted in a 
small jump in holonomy, since the configurations with a 
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T/T, 


0.B7 1. 1.15 1.31 1.51 1.73 1.98 2.27 2.6 



0.B7 1. 1.15 1.31 1.51 1.73 1.98 2.27 2.6 



FIG. 6: Self-consistent value of the holonomy v (upper plot) 
and Polyakov line (lower plot) as a function of action S (lower 
scales), which is related to T/Tc (upper scales). The error bars 
are estimates based on the fluctuations of the numerical data. 


holonomy of 0.5 but with slightly higher density than the 
flat ones, would end up with a smaller free energy. As 
a result we do not get a completely smooth transition 
though that is hidden by the size of the errors as seen in 
Figure and it also means that at S' = 6 the Debye mass 
never goes completely to zero, as shown in Figure and 
the density goes slightly more up also as shown in Fig. 

[i 

When we are in the confined region we observe the free 
energy for a fixed density as a single minimum in the 
middle at = 0.5. As the action S increases, the density 
decrease and it becomes more favorable to have some 
bigger, but lighter dyons, thus shifting the minimum to 
the sides as can be seen in Fig. [T^for S = 6,7, 9. This at 
the same time makes the lighter dyons, more abundant 
than the more heavy dyons. 

Due to the size of the Debye mass, the correlation func¬ 
tions behaves as a liquid with a cutoff at small range. 
We show the case for 5” = 6 in Fig. [D for MM and 
ML. Note the correlation function Cmm vanishes at 
small distances due to the core. The other correlation 
function Cml for ML, displays attraction even at small 
distances, tripling the density at r = 0. The integrated 
number of particles in the region in which the correlation 
function Cml(^) > 1 is 0.50 particles, while for Cmm it 
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FIG. 7: (Golor online). Self-consistent value of the Debye 
Mass M_d as a function of action S (lower scale) which is 
related to T/Tc (upper scale). The error bars are estimates 
based on the fluctuations of the numerical data. Points rep¬ 
resent lattice data from m as a function oiTjT^. 
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FIG. 8: (Golor online). Density n (of an individual kind of 
dyons) as a function of action S (lower scale) which is related 
to TjTc (upper scale) for M dyons(higher line) and L dyons 
(lower line). The error bars are estimates based on the density 
of points and the fluctuations of the numerical data. 


corresponds to 0.34 particles: thus the difference is not 
that large. 


VIII. SUMMARY AND DISCUSSION 

As emphasized in the Introduction, an idea that it 
should be possible to understand confinement (as well 
as chiral symmetry breaking) via statistical mechanics 
in terms of collective coordinates of certain topological 
solitons goes back to 1970’s. Four decades later we now 
have its definite realization in terms of the instanton con¬ 
stituent, the instanton-dyons. 

In particular, by identifying classical interaction be¬ 
tween instanton-dyons uni and including them in direct 
Monte-Carlo simulation of the ensemble, together with 
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FIG. 9: (Color online). Self-consistent value of the Free 
energy density / as a function of action S (lower scale) which 
is related to T/Tc (upper scale) for ro = 2 and A = 1.5. 



FIG. 10: (Color online). (Not-self-consistent in holonomy v) 
Free energy density /, here shown as a function of the value of 
the holonomy (in form of the Polyakov loop P) at S' = 6, 7, 9. 
The lower the action the lower the minimum of the free energy 


one-loop effects in the measure, we calculated the free 
energy as a function of all parameters of the model, such 
as the value of the holonomy, dyon densities, and the De¬ 
bye mass. We then proceed to one-parameter set of its 
minimum, corresponding to dependence on the only left 
variable, the temperature. The results display the de¬ 
confinement transition at a certain density of the dyons. 
The key to this is the volumes of the dyon repulsive cores, 
which scale as an inverse cube of the holonomy. 

One of the key questions usually asked in conjecture 
of such theory is whether the objects we study are suf¬ 
ficiently semiclassical, so let us start our summary with 
answering it. The action per M dyon, Sv, varies in the 
region studied in this work in the range from 2.5 to 3.3. 
Its exponent exp{—Sv) varies between 0.082 and 0.037. 
The input formula we use include classical and one-loop 
effects. By selecting specially tuned A parameter we ba¬ 
sically include the two-loop effects as well. So, we think 
that the accuracy of these expressions is sufficient for our 
purposes. 

Since we perform direct simulations using such expres- 




r 

FIG. 11: (Golor online). Gorrefation function Cij for MM 
and ML for S' = 6 (upper) and S = 9 (lower). In the MM 
case the correlation function vanishes at small distances due 
to the core. 


sions, we have no further approximations, and the ac¬ 
curacy of the results is limited only by statistical errors 
of the Monte Carlo simulations we make. Let us note 
that we did not aimed for high statistical accuracy in 
this work, considering it to be a demonstration of the 
principle. 

What is shown in this work is that the ensemble of 
instanton-dyons, coupled to holonomy, does undergo a 
deconhnement phase transition at certain value of their 
density. It is physically driven by repulsive interactions, 
which enforce “equality” between M and L dyons, broken 
in the dilute regime by their different actions. We see 
how it happens in detail: first by performing multiple 
simulations as a function of all parameters of the model 

- dyon densities, holonomy, the value of the Debye mass 

- and then identifying a co-dimension 1 set of the free 
energy minima, corresponding to physical dyon ensemble 
as a function of the temperature T. All these results 

- the holonomy potential and the mean Polyakov line 
(P(r)), the dyon densities nM{T), nL{T) can and should 
be compared to lattice data. 

This dyon ensemble can be straightforwardly general¬ 
ized to the QCD-like theories with arbitrary number of 
colors and quark flavors. We plan to do larger scale sim¬ 
ulations of those in subsequent publications. 

Finally, let us address a very general question often 
asked: why should one study statistical mechanics of 
some solitons, rather than directly simulate gauge helds 
on the lattice, from the first principles? 
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Quantum field theories have infinitely many degrees 
of freedom, and an understanding of which ones are re¬ 
sponsible for a particular phenomenon is very important. 
Using an analogy to condense matter physics: One can in 
principle do direct simulations of all electrons in a piece 
of metal. And yet, understanding the zone structure, 
location and shape of the Fermi surfaces offer a much 
simpler and more intuitive approaches to metal thermo¬ 
dynamics and kinetics. To a large extent, the same is 
true for quarks in the “zero mode zone” of the topolog¬ 
ical solitons. Now we see that instanton-dyons generate 
confinement as well as chiral symmetry breaking. The 
model we use have only few variables per fmP volume, 
5-6 orders of magnitude less than current lattice simula¬ 
tions. 
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Appendix A: Units and holonomy 


ons, and the corresponding Gross-Pisarski-Yaffe holon¬ 
omy potential [ 2 ] 

VgpyH _ 

3 ^ ’ 

where V 3 is the 3-volume of the box and 

D = 1 — V (-A.3) 

is “dual holonomy” . We proceed in the text to use 
dimensionless units for volume U 3 = T^V^, densities 
P-M = pl = ^, distances rT = x and free en- 

V 3 ^3 

ergy density = /. Potential Vqpy has a minimum 
at trivial holonomy v = 0 and a maximum at confining 
holonomy 1 / = 1 / 2 , thus disfavoring confinement. 

In the next approximation the so called quantum loop 
effects are incorporated. As is well known, they lead to 
a running coupling constant. Thus the action parameter 
(and all others, of course) become a function of the basic 
physical scale given by the temperature T. For example, 
recalling classical instanton action and the asymptotic 
freedom formula 


The main physical quantity of the problem is the tem¬ 
perature T: it defines the magnitude of the = 2ttvT 
(holonomy), the physical size of the dyons and every 
other dimensional parameter of the problem. Yet, pre¬ 
cisely because of its omnipresence in the theory at its 
classical level, dealing only with the dimensionless quan¬ 
tities - e.g. the dyon density normalized as n/T^ - one 
can in zeroth approximation cancel all powers of T. At 
this level, our theory has only dimensionless input pa¬ 
rameters. Most of them - the dimensionless dyon densi¬ 
ties, holonomy and the the Dehye screening mass - will 
be defined selfconsistently, from the minimum of the free 
energy. The remaining input will be the instanton action 
parameter S, used in many plots in the text. 

Standard Euclidean formulation of the gauge theory 
at finite temperature T introduces periodic (Matsubara) 
time T defined on a circle with a period equal to the 
inverse temperature 1/T. The exponential of the gauge 
invariant integral over this circle, known as the Polyakov 
line 

P = ^TrPexp[i /2)dT] (Al) 

which is gauge invariant due to periodicity. Here is 
the 3rd Pauli matrix. 

As a function of temperature its expectation value (P) 
changes from 1 at high T to (near) zero at the decon¬ 
finement temperature Tc- In the simplest SU(2) gauge 
theory we will discuss in this work (P) = cos^ptt), and 
the holonomy parameter (or just holonomy, for short) p 
changes from 0 to 1/2. What remains unknown is the 
physical origin of this potential. 

Perturbatively, the effect of the holonomy is appear¬ 
ance of nonzero masses of quarks and (non-diagonal) glu¬ 


with the power given by the one-loop beta function. If 
so, the semiclassical factors defining the caloron density 
now depend on T, basically as a power 

Pcalorons{T) a (^\ / a rN 

P4 ^\t) ^ 

Since the caloron density has been measured on the lat¬ 
tice at different T, one can test this expression against 
the lattice data. In fact it does work, see Fig.l of Ref. [3], 
which confirms that the topological solitons remain semi¬ 
classical at the temperatures we discuss. 

The next question is the value of the parameter A in 
the expression for S above. Note, that our parameter A 
is proportional to that in multiple other definitions, such 
as Aiattice Or A^g, but is not equal to them. In principle, 
the relation between them is known, and the reader may 
thus ask why we don’t use such relations, obtained from 
the first principles. The answer is pragmatic: we believe 
that the current accuracy of them raised in high power, 
e.g. A^g, is still lower than what was found from the fit 
to the caloron data just mentioned. In other words, the 
measurements of the caloron density is basically the mea¬ 
surements of the high power of A, and they thus provide 
more accurate values than what can possibly be done by 
(much more accurate) measurements but of quantities 
depending on this parameter logarithmically. 

Not surprisingly, in practice the meaning (and the 
value) of A depends on the context in which it is used. 
The fit shown in Fig.l of Ref. [H] corresponds to non¬ 
interacting gas of calorons, and it give 

■^calorons — 0.36Tc (A6) 
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where is the deconfinement transition temperature, 
defined in the same lattice work. If so, the (instanton 
action) parameter is S' « 7.5 at T^. 

In our work we worked out much more sophisticated 
model of the interacting dyon plasma. In this model the 
deconfinement transition happens at a somewhat differ¬ 
ent value of the (instanton action) parameter S « 6 . In 
other words. 


to compensate for the extra 2-kT. The classical interac¬ 
tion between the dyons can at long range be described 
by the same formula for all 

V{r) = ^ f(eie 2 - 2 /tih 2 )-+ TOim 2 -^ (B2) 

9 \ X xj 

X = 2TTvr 


^dyonic plasma — 0.44T). 

This value is assumed in all plots in section [Vll| in which 
our input parameter S is mapped to the temperature T. 

The reader should however keep in mind, that this 
mapping between the input parameter of the model S 
and physical T is provisional, it depends on the model 
itself. No doubt it should be subject to future improve¬ 
ments, both in the lattice data quality used for such fits, 
or the fit itself. In particular, one should include mea¬ 
surements of the instanton-dyons, including the effects of 
their interaction as discussed in the bulk of our paper. 



M 

M 

L 

L 

e 

1 

1 

-1 

-1 

m 

1 

-1 

-1 

1 

h 

1 

1 

-1 

-1 


TABLE I: Quantum numbers of the four different kinds of 
the instanton-dyons of the SU(2) gauge theory. The first two 
rows are electric and magnetic charges, while by h we mean 
the contribution from nonlinear terms including the holonomy 
field. 


Appendix B: Instanton-dyons 

We do not present here extensive introduction on the 
configurations and their history, which can be found in 
literature such as [7]. 

“Higgsing” the SU(2) gauge theory by nonzero VEV of 
A 4 called V leads to two massive and one massless gluons. 
The simplest gauge is the so called regular (hedgehog) 
gauge, in which the color direction of the “Higgs” field is 
at large r along the unit radial vector AJ* —)■ wr™. The 
solutions are 

A 4 = ±fa ^-ucoth(ur)^ 

(r " sinh(ur)) ’ 

where -I- corresponds to the M dyon and — corresponds 
to the M dyon. r is the length in position space. The L 
and L dyon are obtained by a replacement v —>■ 2'kT — v. 
To study the classical interaction of the dyons, a gauge 
transformation is done to make A 4 field point in a spe¬ 
cific direction (normally this is chosen to be A 4 ), which 
introduced a time dependence in the L dyons in order to 


where e, m, h are listed in the Table 1. 

As a result sectors that are completely self-dual or anti- 
self-dual have no interaction, while dyons and antidyons 
of same type attract and dyons and antidyons of different 
type repel. 


Appendix C: The dyon weights in the partition 
function 

The KvBLL caloron partition function [16] has the 
form 


ZkvBLL = 


X 

X 


(27r + i/DTr 12 ) {2TriyTr 12 + 1)^“^ 


5 

3 


(27r:^Tri2 -t- 1 )*^" 1 exp(—V 3 T^^^^z/^:/^)(Cl) 


Taking the limit to very dilute situation we find that all 
powers of Tr 12 not in the exponential cancel, and we end 
with 


ZkvBLL = J d^Z4d^Z2T^c(^^^ (^e ^ 

X (27ri^)T“ (27ri>)“3“ 

A'TT^ 

X exp(-V3r^ —(C2) 


The term in the exponential corresponds to the pertur¬ 
bative holonomy potential. The Diakonov determinant 
which we have included is seen to return to a product of 
the holonomies in the dilute limit 


lim det G = TT 

Tri 2 —>-oo 


(C3) 


By comparison we see that we have to take equation C2 


and divide by equation C3 in order to get the correct 
weight for our partition function. We thus end up with 
the partition function for a M and L dyon given by 


ZrvBLL = 


d^zid^Z2T^C ( — 


8^2 

9 


(27r)®/3 8 ;._ 4 _ 8£_1 

X V 3 ^3 

(47r)2 


(C4) 

















We redefine the constant so the equation is easier to work 
with 
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ZkvBLL = f cf Zi(f Z 2 T^ -7^-^ (] e s' 


(47r)2 
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